Superconducting Fluxon Pumps and Lenses 
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We study stochastic transport of fluxons in superconductors by alternating current (AC) recti- 
fication. Our simulated system provides a fluxon pump, "lens" , or fluxon "rectifier" because the 
applied electrical AC is transformed into a net DC motion of fluxons. Thermal fluctuations and 
the asymmetry of the ratchet channel walls induce this "diode" effect, which can have important 
applications in devices, like SQUID magnetometers, and for fluxon optics, including convex and 
concave fluxon lenses. Certain features are unique to this novel two-dimensional (2D) geometric 
pump, and different from the previously studied ID ratchets. 
PACS numbers: 05.40.-a, 74.60. Ge, 05.60.-k 



Introduction. — A number of authors (see, e.g., H-|l]) 
have recently addressed the longstanding problem of how 
to extract useful work from a fluctuating environment. 
While heat may not be transformed back to mechanical 
work at thermal equilibrium (i.e., without temperature 
gradients (5|), the same restriction does not apply to the 
case of nonequilibrium thermal fluctuations: An asym- 
metric device (like Feynman's ratchet ||^) can indeed rec- 
tify symmetric quasi-equilibrium fluctuations ]^^^ . The 
implications of such a mechanism in transport theory are 
far-reaching: macroscopic currents may arise even in the 
absence of external forces or gradients. 

Recently, a few groups P,0 have made interesting pro- 
posals for quite distinct ways of using ratchets in super- 
conductors. Ref. Q uses Josephson junctions, in SQUIDs 
and arrays, while Ref. M] uses a standard ID potential- 
energy ratchet (e.g., W,M0|) to drive fluxons out of su- 
perconducting samples using simulations like in P,p|. 
Here, we use the new concept of 2D asymmetric channel 
walls to provide an experimentally-realizable way to eas- 
ily move fluxons through channels in AC-driven devices. 
These novel geometric ratchets would allow, for exam- 
ple, the construction of fluxon optics devices, including 
concave/convex fluxon lenses that concentrate/disperse 
fluxons in nanodevices. Further, we also consider the ef- 
fects of collective interactions on ratchets, instead of the 
usual one-particle stochastic transport. 

Simulation. — Computer simulations were performed 
using a new version of the molecular dynamics (MD) code 
used for systems with random ||8| and correlated Q pin- 
ning. We model a transverse 2D slice (in the x-y plane) 
of an infinite zero-field-cooled superconducting slab con- 
taining current-driven 3D rigid vortices that are parallel 
to the sample edge. The samples have very strong, effec- 
tively infinite, pinning except the "zero pinning" central 
saw-tooth- shaped channel (see inset of Fig. 1). In the lat- 
ter, fluxons moved subject to: fluxon-fluxon repulsions 
fyy, an externally- applied AC Lorentz driving force f^, 
forces due to thermal fluctuations fr, and interactions 
with the channel boundaries, treated as infinite potential 



barriers. The thermal force fy was implemented using 
a Box-MuUer zero-mean, Gaussian distributed random 
number generator. We measure all forces in units of /o = 



r2\3 



A , magnetic fields in units of ^o/A , and lengths 



in units of the penetration depth A. Here, $o = hc/2e 
is the flux quantum. Our dimensionless temperature is 
T = /csTactuai/A/o. Thus, away from the channel edge, 
the total force on a fluxon is f = {^ + t + fr = '7V, 
where the total force on vortex i due to other vortices 



F(i) 



^^'^, 



is given by C^ = J2jZi /o Ki{\r, - rj|/A)fy. Here, Ki 
is a modified Bessel function, r^ (v^) is the location (ve- 
locity) of the ith vortex, Ny is the number of vortices, 
fy — {vi — rj)/\ri — Tjl, and we take At = 0.01 and 
r] — 1. Simulations for the same parameters were re- 
peated with different random number seeds, to provide 
for disorder averaging. 

Unless otherwise specified, each figure refers to simu- 
lations conducted in the following way: Initially fluxons 
were randomly placed inside the channel and subjected to 
an alternating current along y — producing a square- wave 
Lorentz driving force along x with F = /l//o = 15 and 
period P = 2r (here, r = 100 MD steps). This means 
that in the absence of thermal noise {T = 0) and channel 
walls, a single fluxon would alternate traveling 15A in one 
direction (e.g., +a;), and then 15A in the other (— .x). The 
period 9A of the horizontal ratchet teeth was such that 
the driving force F = 15 should be sufficient to allow 
for a rectification, or diode effect, because it moves the 
fluxon back and forth through the bottlenecks indicated 
in the inset of Fig. 1. Ten runs, of 250000 MD steps each, 
were used to find the average current and standard devi- 
ation for each plotted point. When a standard deviation 
is not shown, it was significantly smaller than the plot- 
ted point. All samples had a 2D geometric ratchet made 
of asymmetric walls. The simulated sample was an 18A 
by 18A square, with periodic boundary conditions. The 
channel is 7A wide, with saw-teeth of period 9A (four 
teeth per unit cell) and slope 1/3 (except in Fig. 5, and 
in other runs to be discussed elsewhere, where geometry 
was varied). This leads to a bottleneck that is one A 



wide. This construction is a novel 2D ratchet that uses 
geometry, which governs the system potential energy, to 
cause rectification. 

Temperature dependence of the fluxon rectification. — 
Figure 1 shows the rectified fluxon velocity {v) = 
J2i " "^i/^v versus temperature. The dimensionless tem- 
perature T was varied over two orders of magnitude, from 
1 to 100. These simulations clearly indicate an opti- 
mal or "resonant" temperature regime in which the DC 
fluxon velocity is maximized by the fluxon pump or diode. 
This optimal temperature regime can be explained as a 
trade-off between allowing the fluxons to fully explore 
the ratchet geometry (i.e., T must not be too low) and 
washing out the driving force (i.e., T must not be too 
large). At low temperatures, the alternating driving force 
will cause a fluxon to migrate to the center of the chan- 
nel and no longer be impeded or assisted by interactions 
with the geometry; while at high temperatures the driv- 
ing force becomes irrelevant and thus the fluxon is no 
longer pushed through bottlenecks regularly. The (w)(T) 
curve for many randomly placed fluxons is similar to the 
single fluxon case, but the magnitude of {v){T) decreases 
when the fluxon density B increases. This is because 
the repulsive force produced by a large number of flux- 
ons act to restrict each other's motion. Note that at the 
optimal temperature, the fluxons travel nearly 15 A ev- 
ery 100 MD steps. The general shape of the curve was 
established with two fine scans with one and fifty fiux- 
ons. Conformity of other field strengths to this curve was 
judged by performing a coarser scan, then connecting the 
points with a spline. The small dip in the (v) for low T 
will be discussed separately (see Fig. 3 below). 

Driving period dependence of the flux pump. — Fig. 
2(a) shows {v){P) for six different combinations of driv- 
ing force amplitudes F and temperatures T for a single 
fiuxon. Interestingly, as the period P of the driving force 
was varied from low to high (from P = 20 MD steps 
to 2000 MD steps) no optimal peak in (v) was discov- 
ered. This is in contrast with previous work on ratch- 
ets, which provide a peak in (v) versus either P or fre- 
quency. Instead, we find that while a high frequency driv- 
ing force yields a very small ratchet velocity, the velocity 
quickly converges to a nearly stable value with increas- 
ing P. This result corresponds to the idea that a fiuxon 
must be forced through the bottleneck of the geometric 
2D ratchet in order for it to be a rectifier. At high fre- 
quencies, a fluxon docs not travel far enough to interact 
with the bottleneck. Once the frequency is low enough to 
force it through a bottleneck, however, the net velocity 
changes little by forcing it through further bottlenecks. 

All of the simulations in Fig. 2(a) were conducted 
at relatively low temperatures to allow the effects of 
the driving force to dominate. In addition to trying 
F ~ /l//o ~ 15, the value typically used in this pa- 
per, two combinations of driving force and temperature 



magnitudes were tried that both also had an F/T ratio 
of 3/1. While varying P at the low-T regime had little 
effect upon the period dependence curve, varying F very 
clearly did have an effect: The higher F, the greater {v). 

Fig. 2(b) shows the dependence of (v) on P, for T — 
F/3 = 5 with many fiuxons instead of one. Initially, fifty 
randomly placed fluxons were simulated. As with the 
single fluxon case, a high frequency driving force resulted 
in a low (v) because the fluxons were not being forced 
through the ratchet bottleneck. Increasing P produced 
a rapidly plateauing (v). Once the shape of the curve 
was determined, similar sets of observations were made 
of different fluxon densities B at six different periods. 
These simulations of many fluxons demonstrate that, as 
with one fluxon, there was no optimal peak. 

Fluxon pump effect versus field. — By varying the den- 
sity of fluxons B, we found a maximum in {v) at very low 
temperatures (see Fig. 3). This is due to an enhanced in- 
teraction with the ratchet geometry: at extremely low T 
the driven fluxons eventually work their way to the mid- 
dle of the channel and then no longer interact with the 
ratchet. A small increase in B can drive fluxons from 
the middle of the channel, causing these fluxons to in- 
crease their interaction with the sawteeth, without sig- 
nificantly increasing the resistance due to fiuxon-fiuxon 
interactions. Aside from this low-T resonance, increas- 
ing B slowly decreases (u), because, as the fluxons repel 
one another, at high densities they block the bottleneck. 
At high temperatures (e.g., T = 50) the ratchet velocity 
decreases slowly with increased fluxons. In this case, the 
slight maximum in (v) found at low T is not present be- 
cause the resistance due to fiuxon-fiuxon interactions is 
always larger than the gains from increased exploration 
of the geometry of the ratchet. 

Flux pump response to the driving force. — Increasing 
the amplitude F of the square-wave driving force in- 
creases the average velocity (w) of a fiuxon in the ratchet, 
as shown in Fig. 4. This monotonic increase in (v) ver- 
sus F for 2D geometric ratchets is different from the peak 
observed in ID potential ratchets (e.g., see Fig. 2 of [Q). 
This makes it fairly easy for any particular 2D geometric 
ratchet to be continuously tuned to desired values of (w) . 
This would be an advantage of 2D geometrical over ID 
potential ratchets for potential practical designs of de- 
vices. A reason for their different response is because the 
resistance to motion due to successive bottlenecks does 
not add like the resistance due to successive potential 
barriers. If a fluxon happens to remain for some time 
in the center of a bottleneck, unlikely to happen at high 
enough temperatures, then it will not be rectified; while 
the potential barrier rectifies at every period. Increased 
fluxon density again increases resistance, reducing (v). 

Geometry dependence. — To determine the generality of 
the simulation results presented here with respect to var- 
ied ratchet geometries, in Fig. 5 we compare (w)(T) for 
fifty fiuxons in two alternate geometries with the (w)(T) 



shown in Fig. 1 (open squares in both figures) . Increasing 
the slope, and correspondingly decreasing the period, of 
the saw-teeth clearly shifts the optimal temperature to 
lower values. Despite this shift, however, the qualita- 
tive behavior of the first geometry studied carried over 
to other ratchet slopes and periods. Additional samples 
and parameters, including wider bottlenecks, were also 
studied (but not discussed here due to space limitations), 
giving consistent results. 

Discussion. — Our 2D geometric ratchet has some 
properties which are similar to previously studied ID 
potential ratchets. Like the simpler ID cases, the 2D 
ratchet showed a "resonance region" with a maximum 
in (v), and (outside of a low-T anomaly) a decreased 
(v) with increasing B. Unlike the ID cases, however, 
the ratchet velocity plateaus with decreased driving fre- 
quency, instead of displaying a peak. Also, 2D geomet- 
ric ratchets are easier to implement experimentally since 
these do not require carving very many inclined asym- 
metric grooves, as with a ID potential ratchet |7|. More- 
over, as pointed out above (fig. 4), {v) can be continu- 
ously tuned by varying F, for the 2D geometric ratchets. 

Concluding remarks. — Superconducting devices using 
stochastic transport were studied because they provide 
a convenient way to move fluxons through channels. At 
non-zero temperatures, the asymmetric geometry of the 
ratchet saw-teeth automatically converts applied AC in- 
puts into a net DC motion of fluxons. Thus, the simu- 
lated device serves as an excellent fluxon rectifier. Our 
device could be easily made experimentally either by elec- 
tron beam lithography, irradiation, evaporating layers, 
or by chemically etching a channel. The central channel 
would have very weak pinning, while the rest very strong 
pinning. By coupling two ratchets that rectify in oppo- 
site directions (see inset of Fig. 2(b) and Ijl^), fluxon 
lenses that could either (c) disperse or (d) concentrate 
fluxons in chosen regions of a sample can be created. 
More complex fluxon optics microdevices could be built 
similarly. Corners (e) can also be constructed. Such re- 
markable devices, and modifications of them, would al- 
low the transport of fluxons along complicated nanofabri- 
cated channels: a microscopic network of fluxon channels 
in superconducting devices. This could be very useful to 
get rid of unwanted, trapped flux in SQUID magnetome- 
ters, and also to move fluxons along channels in devices 
pi. These promising concepts are largely unexplored, 
and constitute an open and potentially useful area. 
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FIG. 1. Rectified average fluxon velocity {v), which can be 
measured as a voltage, versus temperature T for the ratchet 
geometry in the inset. The magnetic field B is directed out of 
the figure. J = Jy is a vertically applied alternating current 
square-wave that drives the fluxons back and forth horizon- 
tally along the channel. The number of fluxons is: 1 (0)i 25 
(O), 50 (□), 75 (A), 100 (x), 150 (+), 250 (*). 



FIG. 2. Average fluxon velocity (v), versus driving force 
period P for (a) one fluxon and (b) many fluxons. For very 
low periods little rectification occurs, and for higher periods, 
(v) slowly increases. The four F — 15 curves in the middle of 
(a) show similar behavior and correspond to (top to bottom 
at r = 10): T = 1 (•), T = 2 (O), T = 3 (x), T = 2 (+). In 
(b), similar (w) occur for 50 (O), 100 (□), 150 (A), and 250 
(x) fluxons. Inset: Schematic diagram for concave (c) and 
convex (d) fluxon lenses that disperse/focus fluxons from/at 
their centers, (e) shows a corner unit. 



FIG. 3. (v) versus the number of fluxons, for several values 
of the temperature T. 



FIG. 4. (v) versus magnitude F of the driving force. The 
number of fluxons for the low-T (T — 5) cases are: 1 (0)i 
50 (O), and 100 (□). 50 fluxons were used for the high-T 
(T = 50) curve (A). 



FIG. 5. (v) versus T for different ratchet geometries, char- 
acterized by {S,Pst), where S is the slope with respect to x, 
and Pst the spatial period of the sawtheeth. The curves shown 
correpond to (1/3,9) (D), (2/3,9/2) (A), and (1,3) (Q)- 
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